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ABSTRACT 

We derive galaxy source counts at 70 and 160 /xm using the Multiband Imaging Photometer for 
Spitzer (MIPS) to map the Chandra Deep Field South (CDFS) and other fields. At 70 fim, our 
observations extend upward about 2 orders of magnitude in flux density from a threshold of 15 mJy, 
and at 160 /im they extend about an order of magnitude upward from 50 mJy. The counts are 
consistent with previous observations on the bright end. Significant evolution is detected at the faint 
end of the counts in both bands, by factors of 2-3 over no-evolution models. This evolution agrees 
well with models that indicate most of the faint galaxies lie at redshifts between 0.7 and 0.9. The new 
Spitzer data already resolve about 23% of the cosmic far-infrared background at 70 //m and about 7% 
at 160 /im. 
Subject headings: infrared: galaxies - galaxies: evolution - galaxies: statistics 



1. INTRODUCTION 

The cosmic infrared background (CIB), relic emission 
of the formation and evolution of the galaxies, peaks in 
the far-infrar ed (FIR) i n the 60-200 um wavelengt h range 
llPuE jet et all r i996: Ha i|ser et all [19 98: Lagach eet all 
11999"; "Gispert et al. 2000; H auser fc D wck 2001). In the 
local universe only about a third of the cxtragalactic 
emiss ion is released in the FIR (jSoifer fc Neugebaueil 
Il991[) . However, the CIB FIR peak accounts for more 
than half of the total optical/infrared background, indi- 
cating strong evolution of galaxy properties toward high 
FIR output in the past. Characterizing the galaxies re- 
sponsible for most of the CIB is therefore an important 
goal of cosmological surveys. Galaxy counts (or number 
counts) provide a powerful tool to investigate the evolu- 
tion of the galaxies and their contribution to the CIB. 

The cryogenic infrared space missions IRAS (In- 
frared Astronomical Satellite) and ISO (In f rared 
Space Observatory) (" Genzel fc Cesarskvl 1200(1: iDolel 
l2003[ for reviews) provide d important data on source 
counts at 60 ixui ( Hacking & Soifer "1991J) and at 
100 /im IjRowan-R obinson et al. 1 986) , and more 
recentl y at 90 um JK awara eLaUll^8; Efsta tiiiou et al.l 
J .2000: iLinden-Vornle et alJ 
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Kawara et al." '2004") and 170 /im ("Kawara et al.l 
Pyge t et al. 1999; Matsu hara et al. 2000; Dol e et all 
200H iKawara et all 1200^ . Mid -infrared (MIR) o bser- 
vation with ISOCAM at 15/xm l|Elbaz et al.lll999() are 
also of great interest since they are believed to resolve a 
significant fraction of the CIB into sources IjElbaz et al.l 
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I2002|) . The Spitzer Space Telescope l)Werner et alJl2004|) 
provides the ability for much deeper and wider-area 
surveys from 3.6 to 160 /im. This paper investigates 
source counts at 70 and 160 /^m from Spitzer. A 
compan ion paper addresses t he MIR source counts at 
24 /im IjPap ovich et al.ll2004j) . The three-band source 
counts are th e basis of new phenomenological models by 
lLagacheet"al.ll|20d^ . 



2. OBSERVATIONS AND DATA REDUCTION 

Observations were carried out with the Multiband 
Imaging Photometer for Spitzer ( Ricke ct al. 2004) in the 
Chandra Deep Field South (CDFS) and the Bootes field 
corresponding to the NOAO D eep Wide Field Survey 
NDWFS (Jann uzifc Devlll999|) : We also used an engi- 
neering MIPS observation of the Marano field. The ob- 
servational mode (scan map) provides multiple sightings 
of each source, typically 10 and 60 at 70/im in the Bootes 
and CDFS, respectively. However, at 160 /xm, the num- 
ber of sightings is onl y typically 2 in Boote s and 12 in 
CDFS. See TablelTland rPaDOvich et all (J2^ for details. 

The data were re duced with the Data Analysis Tool 
IjGordon et alJ l2004'). from the raw data (ramps) to the 
final coadded mosaics. The illumination corrections were 
derived from the data themselves. At 70 /xm, the data 
have been median-filtered in the time domain before mo- 
saicking. Note that data from Ge:Ga detectors are al- 
ways challenging to process; but with MIPS, most of 
the difficulties are overcome with frequent calibrations 
(stimulator flashes), that track the responsivity varia- 
tions. Nevertheless, the noise properties at faint fluxes 
are still being investigated at both 70 and 160 fim. In 
this work, we will adopt conservative detection limits. A 
future paper will address extracting the ultimate sensi- 
tivity from these data. Sample images in the CDFS are 
shown in Figure ^ 

3. PHOTOMETRY & CATALOGS 

To control the sample and the selection function, we 
accepted source detections only where the redundancy 
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TABLE 1 
Log of Observations 



Field 


MIPS 




70 Atm 






160 fim 




Name 


AOT'' 


Area 


^int 


Scut'' 


Area 


tint" 


Scut'' 


a 




[deg'] 


[S] 


[mjy] 


[deg'] 


[s] 


[mJy] 


Bootes 


mod. 


8.75 


40 


80 


7.70 


8 


... <^ 


Marano 


slow 


0.42 


100 


25 


0.31 


20 


50 


CDFS 


slow 


0.67 


600 


15 


0.54 


120 


50 



^ see Papovich et al. (2004) for details on fields 

^ scan map mode 

'^ per sky pixel 

'^ flux density at which catalog was cut 

" field data not used at 160 fira 



was high (typicaUy 80% or more of the mean weight), 
avoiding the edges and the noisiest areas of the images. 
The r esulting positions were fed to DAOPHOT IjStetsonI 
[19871) in IRAPS for PSF fitting. We checked that the 
residual maps were indeed free of sources. 

At 70 /xm the photometric calibration is derived from 
many observation campaigns, and its uncertainty is con- 
servatively estimated at the order of 20%. We use only 
detections at 15 mJy and brighter in the CDFS, 25 mJy 
and brighter in Marano, and 80 mJy and brighter in the 
Bootes field. These flux levels are determined using the 
sharp decrease in the counts due to the incompleteness 
effect. At 160 /im, the calibration is based on a combi- 
nation of observations of standard stars, asteroids, and 
comparisons with measurements with other FIR missions 
{ISO, CO BE, and modeling including IRAS measure- 
ments). It is also estimated to be accurate to about 
20%. We have included in our counts only objects of 
50m Jy and brighter in the CDFS, and Marano (levels 
determined using incompleteness as at 70 /im). Because 
of the low redundancy level of the 160 /im data in the 
Bootes field, we postpone using it for a later paper. 

Catalogs were produced separately at each wavelength; 
surveys at each wavelength are thus unbiased. Detected 
FIR sources in the CDFS sort as follows. For sources 
selected at 70 /im: 92% have a 24 /im ID; 54% have a 
160 /im ID (same with 24 and 160). For sources selected 
at 160 /Ltm: 98% have a 24 ^m ID; 43% have a 70 /im ID 
(same with 24 and 70). 

4. SOURCE COUNTS 

At 70 /Lim, 131 sources were detected down to 15 mJy 
in the CDFS, 55 sources down to 25 mJy in Marano, and 
117 down to 80 mJy in the Bootes. At 160 /im, down to 
50 mJy, 123 sources were detected in the CDFS, and 89 
sources in Marano. The source density corresponds to 
about 150 beams per s ource at 70 /xm and 15 at 160 /im, 
using the definition of iHelou fc BeichmanI l)1990f) . 

Source counts are given in integral form (Figure \^ 
and|3^) and differential form, divided by the Euclidean 
component (Figurel^b and^D) at 70 and 160 /im, respec- 
tively. Notice that we did not correct for incompleteness. 
Error bars on counts are \a Poisson uncertainty. Bins 

* IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under the cooperative agreement 
with the National Science Foundation 



with less than four sources have not been displayed for 
clarity, since their significance is low, with uncertainties 
of 50% or higher. Also, photometric uncertainty have 
been displayed at high flux only. In order to visualize 
the contribution from each field, source counts have not 
been merged and have been overplotted. One should 
keep in mind that MIPS source counts will eventually go 
deeper and will be corrected for incompleteness. 

The observed fields nicely complement each other in 
terms of area and depth. This allows us to probe a flux 
range of almost 2 orders of magnitude at 70 /im. One 
order of magnitude is covered at 160 /xm. It is possible 
to check consistency and the cosmic variance in the com- 
mon flux density range. At 70 /im, in the range 25 to 
100 mJy where three fields overlap in flux density, the 
differential counts are almost consistent within the er- 
ror bars. At 160 /im, in the range 100 to 300 mJy, the 
differential counts are consistent within the error bars. 
At both wavelengths, number counts in CDFS appear 
consistently lower than in Marano. 

5. DISCUSSION 
5.1. 70 ^im 

The MIPS 70 /im counts show a gre at cons istency with 
the IRAS 60 /im counts of iLonsdal e et al.l (^990) con- 
verted at 70 /im using vq^Sqq — v-jqS^q. 

A selection of recent models is shown in Figure |2 in- 
cluding a non-evolution scenario. The most striking re- 
sult is the strong excess of MIPS 70 /im sources compared 
to the non-evolution model: a factor of 3 at around 20- 
30 mJy. Strong evolution had been reported previously 
at 60 and 90 /im, and these data provide unambiguous 
confirmation. 

Two models lie close to the data: 
iKing fc Rowan-RobinsonI l|2003D and iLagache et alJ 
l|2003[ 120041) . These models, developed to fit observables 
mostly from IRAS, ISO and SCUBA surveys as well 
as the CIB spectral energy distribution (SED), are 
based on a strong evolution of luminous (and ultra- 
luminous) infrared galaxies (LIRG and ULIRG). The 
latter model predicts a peak in the redshift distributions 
of re solved sources at 70 /im near z ~ 0.7 fDolc ct aO 
120031) . Figure 0^ shows the galaxy contribution to the 
differential counts, a s a function of redshift, from the 
ILagache et al.l l|2004l) model. Between ~5 and --100 
mJy, sources at 0.7 < z < 0.9 contribute the most to the 
counts. At brighter fluxes (reached by IRAS and ISO), 
contributions from local galaxies are more important. 

The source counts integrated at 70 /xm correspond to a 
brightness of 0.022 MJy/sr or 0.95 nW/m^sr. The value 
of the CIB at this wavelength is not known accurately 
owing to contamination by zodiacal light. If we us e the 
CIB value from the model of ILagache et alJ |2004), the 
MIPS counts show that at 70 /xm about 23% of the 70 /im 
CIB is already resolved. 

5.2. 160 ^im 

At 160 /im, t he comparison w ith the ISO FIRBACK 
170 /im survey l|Dole et al.ll200lD shows that the counts 
are consistent and within the error bars in the whole 
common range, 180-300 mJy. Other ISO ob servations 
(jMatsuhara et alJ 120001: iKawara et all l2004|) agree as 
well. 
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Figure 121 overplots the same models as at 70 /im. The 
evolution detected at 170 /im is confirmed at 160 ^m, 
down to fainter levels. At about 100 mJy, an excess 
of sources by more than a factor of 2 is observed com- 
pared to a non-evolution scenario. Interestingly, as for 
70 u m, the observed evolut ion i s better fitted by the mod- 
el s oflLaeache et aLl 1)2004)) and lKing fc Rowan-Robinsonl 
foes'). The observed slope also agrees with iMouldl 
003). 

We have constrained the bright end of galaxy num- 
ber counts at 170 /im by using data from the ISOPHOT 
Serendipity Survey (ISOSS). ISOSS provides a total sky 
coverage of 15% and is virtually complete at a flux den- 
sity level of 5*170 /im = 50 Jy. Based o n all opticall; 
identified galaxies detected by ISOSS ('Krause ' '2002_ 
iStickel etal. 2004) we have derived an integral number 
density of n(5'i 70 ^m > 50 Jy) — 14±3gal.sr~^ galaxies at 
high Galactic la titudes. This p oint is perfectly matched 
by the model of iLagache et all ([2003, 2004). 

The most striking result of the models that fit the num- 
ber counts is the existence of two regimes in flux den- 
sity. In the ISO range (fluxes above 200 mJy) most of 
the sources contributing t o the counts are l ocal: this is 
confirmed by observation IjPatris et alJl2003(l . At fainter 
fluxes, between ~10 and 200 mJy, the counts should be 
dominated by a population located at redshifts between 
0.7 and 0.9. 

The source counts integrated at 160 /xm correspond 
to a brightness of 0.07 MJy/sr or 1.4 nW/m^/sr^_rhe 
GIB value at this wavelength is 1 MJy/sr iLagache et al.l 
|2000); the MIPS counts show that at 160 fim about 7% 
of the GIB is resolved. Since these counts are prelimi- 



nary, are not corrected for incompleteness, and are sub- 
ject to cosmic variance, we anticipate that the actual 
value might be higher. 

5.3. Concluding Remarks 

The first MIPS far-infrared source counts, spanning 
about 2 orders of magnitude in flux density at 70 /im 
(and one at 160 /xm), are consistent on the bright end 
with previous observations, and show on the faint end 
unambiguous evolution. Models predict that most of 
the sources lie at z ~ 0.7 with a t ail up to z '^ 2 
IjDole et al.ll200l ILagache et al.ll2003ll20o|) . This work 
and companion paper s about source counts at 24 ^m 
IjPapov ich et al.' '2004*1 . about confusion at 24, 70 and 
160 fim (Dole ct al. 2004), and about the interpretation 
of these new data fro m the Spitzer cosmological surveys 
IjLagache et al.ll2004|) . bring new light on the statistical 
properties of galaxies in an unexplored regime in flux 
density, and likely in a critical region of redshift space 
(up to redshifts z - 2) in the FIR (,Egami ct al...200l 
ILeFloc'het7ani200l . 
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Fig. 1. — MIPS observations of the Chandra Deep Field South at 70 and 160 fim. The field covers an area of 25' X 1°. 
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Fig. 2. — Source Counts at 70 /im with no correction for incoin'pl^j^ctnys.y Red Star: CDFS; Blue Diamond: Marano; Black Square: 
Bootes Field; Black Diamond: IRAS 60 fiva counts from Lonsdale et al. (1990) converted at 70 ^m. Top panel (a): Integral Source 
Counts. For clarity, photometric uncertainty is only shown for Sro > 100 mjy. Bottom panel (b): Differential Source Counts. Models are 
also plotted: long d ash: Devriendt & Guiderdoni (2000); dash-dot: Charv & Elbaz (2001); dash: IBalland et alJ l|200,'^1 : dash-dot-dot-dot: 
IKing &z Rowan-Robinson. L2003il : solid line: .Lagache et al.. L2004) : plus: No Evolution Model. 
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Fig. 3. — Source counts at 160 /im with no cor r ection for incompletSiesq.rffifei]Star: CDFS; Blue Diamond: Marano field; Black Diamond: 
ISO FIRBACK 170 /xm counts from lDole et al.l 120011) . Top panel (a), integral Source counts. For clarity, photometric uncertainty is only 
shown for 5160 > 200 mjy. Lower panel (b): Difi^erential Source counts. Models are also plotted with the same symbols as Figure 1^ 
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Fig. 4. — Cumulative contributio ns to the Differential Source Counts at 70 fim and 160 /^m of galaxies, as a function of redshift, form 
the model of ILagache et all 120041) . Top panel (a): 70 ^m. Bottom panel (b): 160 fim. Symbols are described in Figures 1^ and 1^ 
Galaxies contributing to the counts at Redshifts 0.1 (dot), 0.3 (dash), 0.7 (dot-dash), 0.9 (dot-dot-dot-dash) and 1.1 (long dash). The total 
contribution is the upper solid line. 



